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Universality in invariant random-matrix models: Existence near the soft edge 
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We consider two non-Gaussian ensembles of large Hermitian random matrices with strong level 
confinement and show that near the soft edge of the spectrum both scaled density of states and 
eigenvalue correlations follow so-called Airy laws inherent in Gaussian unitary ensemble. This 
suggests that the invariant one-matrix models should display universal eigenvalue correlations in 
the soft-edge scaling limit. 
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I. INTRODUCTION 


Unitary-invariant random matrix models appear in 
many physical theories including nuclear physics, string 
theory, quantum chaos, and mesoscopic physics. They 
are completely defined by the joint distribution function 


P[H] = ^exp{-TrU[H]} (1) 

of the entries of the N x N Hermitian matrix H. In Eq. 
(1) the function V [H] is referred to as “confinement po¬ 
tential” , and it should provide existence of the partition 
function Z]\f. Remarkable feature of this random ma¬ 
trix model is that, under certain conditions, a particular 
form of the confinement potential exerts no influence on 
the local eigenvalue correlations in the bulk scaling limit. 
More precisely, there is a class of strong even confining 
potentials V (e), increasing at least as fast as |e| at in¬ 
finity, for which the two-point kernel in the bulk of the 
eigenvalue spectrum follows the sine form in the large- 
limit IHl: 


A^bulk ) 


sin [tt (s — s')] 
TT (s — s') 


( 2 ) 


This striking property, known as local universality, 
leads to the conclusion about universality of arbi¬ 
trary n-point correlation functions i?„ (si,..., s„) = 

det [K (si, Sj)] ■ ^ (n > 1) in the local regime. In con¬ 

trast, the global characteristics of eigenspectrum, like 
density of states or one-point Green’s function, display 
a great sensitivity to the details of confinement potential 

I- 

Less is known about eigenvalue correlations near the 
soft edge which is of special interest in the matrix models 
of 2D quantum gravity ^. In the early study [|j the be¬ 
havior of the density of states near the tail of eigenvalue 
support has been explored. Authors of Ref. Q showed 
that there is a universal crossover from a non-zero den¬ 
sity of states to a vanishing density of states which is 
independent of confining potential in the soft-edge scal¬ 
ing limit. Whereas the universal behavior of the density 


of states in the soft-edge scaling limit has been proven, 
the (supposed) universality of n-point correlations was 
not yet considered. 

The problem we address in this work is: Whether the 
eigenvalue correlations in ensembles of large random ma¬ 
trices also possess a universal behavior in the soft-edge 
scaling limit? To provide an answer to this question we 
first quote some results for Gaussian unitary ensemble 
(GUE) that has received the most study near the soft 
edge, and then we turn to the consideration of eigenspec- 
tra of two strongly non-Gaussian ensembles of random 
matrices associated with quartic and sextic confining po¬ 
tentials. 

In the soft-edge scaling limit, GUE is characterized by 
the Airy two-point kernel [Q 

Ai (s) Ai'(s') - Ai (s') Ai'(s) 
A.gue(s,s) =-;-, (3) 

c — o' 


whose spectral properties have got the detailed study in 
Ref. 0. As a consequence of Eq. (3), the scaled density 
of states, unlike in the case of bulk scaling limit, cannot 
already be taken as being approximately constant, and 
changes in accordance with Airy law: 


i^GUE (s) = ( ^ Ai (s)) -s[Ai(s)]^ (4a) 


with asymptotes 


i^GUE (s) = 


\/H cos(4|s|’^'’^/3) 


47r|fi 


Stt 


ijexp (-4 s 3/^/3) , s^-boo. 


(4b) 


Our following treatment of non-Gaussian random ma¬ 
trix ensembles with strong level confinement will be built 
upon the orthogonal polynomial technique ^ allowing 
to express the two-point kernel for the random-matrix 
ensemble defined by Eq. (1) through the polynomials 
Pn {s) orthogonal on the whole real axis with respect to 
the weight exp{—2U (e)}. We fix the polynomials Pn 
satisfying the three-term recurrence formula 


^Pn — ^n-\-lPn-\-l “b ^nPn — 1 


(5) 
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to be orthonormal, 


r*+oo 


dePn (e) Pm (e) exp {-2V (e)} = 6n 


( 6 ) 


Thereafter we shall be interested in the behavior of the 
wavefunction tpn near the soft band edge Dn in the limit 
n = N ^ 1. In this case the endpoint of the spectrum 
Dpf = 2aM, where 


Under these conditions the two-point kernel reads as 


_ , ,, i’N {s') i^N-i (e) - tpN (e) t/’w-i (p) 

Kn (£,£ ) = aN - ) - 


(7) 


where 0 qn = kN-i/kN [kN is a leading coefficient 
of the orthogonal polynomial Pn (£)], and the “wave- 
functions” ipN {£) = Pat (£) exp {—U (e)} have been in¬ 
troduced. Inasmuch as our concern is with the matrices 
of large dimensions, N ^ 1, the only asymptotics of the 
“wavefunctions” ipN are needed, and also a meaningful 
scaling limit should be constructed. Quite generally, this 
can be done by passing from initial energy variable e to 
a new scaled variable s that remains finite as N oo: 
e = e{N,s) = Eg- Then the scaled two-point kernel is 
determined by the formula 

K{s,s')= lim Kn {es,£s') (8) 

N^oc as 


II. QUARTIC CONFINEMENT POTENTIAL 


We choose the quartic confinement potential in the 
form V (e) = In this case the differential equation 
for ipn (e) [index n is arbitrary positive integer] can be 
obtained by the Shohat’s method 11 1^: 


de^ 


Ipn (£) - 


de 


Xmpn (e) 


de 


(e) + Qn {£) V'n (^) = 0, 
(9a) 


“"^=( 12 ) + ( 11 ) 

and, 

^N {£) = 2a% + + O (12) 


Let us move the spectrum origin to its endpoint Dn, 
making replacement e — Dn + t, and denote ipN (t) = 
tpN {£ — Dn) ■ It is straightforward to show that this 
function obeys equation 

(t) - 18D%t ■ ipN (t) = 0 (13) 


in the asymptotic limit N ^ 1. When deriving 

we supposed the characteristic energy scale U (N) = 
^ -1 ^ 


dln^N {t) /dt 


of the variation oi ip n (Q to be much 


smaller than the band edge Dn- 

Solution to Eq. (13) can be written through the Airy 
function y (x) = Ai (x) satisfying the differential equation 
y" {x) - xy (x) = 0: 


^N{t) = XNA-X{t-{l8D%f^y . (14) 


One can check that the condition ty (N) ^ 

O <C Dn is fulfilled. The coefficient Xn enter¬ 

ing Eq. (14) still remains unknown. 

To compute the two-point kernel, Eq. (7), we have 
to correctly determine the asymptotic behavior of the 
ipN-i (t). This can be done by means of the asymptotic 
analysis of the exact relation Eq. (10), which in the large- 
N limit comes down to 


^ \ 2 I 2 I 2 

ipn (e) — + a„ + e , 


(9b) 


1pN-l (t) = V'AT (t) + 


1 d 
3D^ dt 


IpN (t) . 


(15) 


Qn (,£) 



4e^ \ 

>fn {£)/ 


+ 4a 


2 

n 


X 


(e) T’n-I (e) -b 1 - 4 a^e^ - 4 e^ 



(9c) 


It is convenient to define the soft-edge scaling limit as 




Dn + 



(16) 


Then the two-point kernel, Eq. (7), and the density of 
states K {£s,£s), are given by the formulas 


Here a„ is the recursion coefficient entering the corre¬ 
sponding three-term recurrence formula for the given set 
of orthogonal polynomials. Also, the following exact re¬ 
lation takes place: 


K (Es, Eg') — X'n 



1/3 

ATgue (s, s') , 


and 


1pn-l (e) = 


i’n (g) + i’n (e) [U' (g) + 4ga^] 
4a„(/>„ (e) 


( 10 ) 


V {Eg) = X% 



1/3 

J^GUE (s) , 


(17a) 


(17b) 
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respectively. The latter expression provides a possibil¬ 
ity to determine the unknown constant Xn by htting the 
soft-edge density of states, Eq. (17b), to the bulk density 
of states Q 


/ N _ 

^ulk (‘^s) — ' 

TT 


\Dn 


l-k2 


(— 

\Dn 


(18) 


taken near the endpoint of the spectrum, Eq. (16), pro¬ 
vided 1 s Equations (18), (17b), (16), and 

(4b) yield the value X% = . Now, making use 

of the Eqs. (17a), (16), and (8), we arrive at the follow¬ 
ing expression for the two-point kernel in the soft-edge 
scaling limit: 


Ksoh (s, s') = Kque ( s , s') . 


(19) 


Thus, we conclude that the two-point kernel and the 
density of states, computed for the random matrix en¬ 
semble with quartic confining potential in the soft-edge 
scaling limit, coincide exactly with those for GUE. 


III. SEXTIC CONFINEMENT POTENTIAL 

Now we turn to another ensemble of random matri¬ 
ces which is characterized by the confinement potential 
V (e) = Corresponding wavefunctions tpn satisfy 

the same differential equation, Eq. (9a), but with Q 




— In ipr, (e) 


+alipn (e) (fin-i (e) - ( e® + 7r„ (e) - — ) 7r„ (e) 


For n = N ^ I the recursion coefficient |n 
1/6 


UN = 


((1) [1 + 0 («-»)], (22) 


and 


ipN (e) = Qa% + {e^ + 2a%) + O (^N , (23) 


ttn (e) = a%£ (e^ -I- 3a%) + O (^N . 


(24) 


Introducing the shifted energy variable, e = Dn +1, we 
are able to rewrite the differential equation (9a) for the 
function 'ipN (t) = if )n {e — Dn) in the form 


^ poc: ^ 

{t)-—Dlt.^l^N{t) = Q. 


(25) 


assuming that the characteristic energy scale ty (N) = 
^ -1 ^ 
dhiipN it) /dt of the variation of 'ipN (t) = 

'ipN {e — Dn) is much smaller than the band edge Dn- 
Solution of Eq. (25) takes the form 


A. ((^)‘'%), 


( 26 ) 


with coefficient Xn that will be determined later by the 
same fitting arguments. The assumption ty (TV) Dn 
is obviously fulfilled. 

To get the asymptotic behavior of tpN-i (t) in the 
large-A limit, we simplify Eq. (21) to 

16 d 

^N-l (t) = ^N {t) + —$N (t) . (27) 


(2^2 + , 


' Vn (e) 


®n-|-l) ) 


(20a) 


and 


7r„ (e) = ale {al_i +al + -k e^) , (20b) 


'-Pn (e) — a„+i (^n+2 ffi ®rt+l + ®n) 


+ On (On +1 + On + 0 ^_i) + (o^+i + O^ -|- e^) . (20c) 

Also, the following relationship holds for arbitrary n: 

, '^I(e) + A (e) [V' (e) + 7r„ (e)] 

Vn-l (EJ = -y--. (2ij 

a„(p„ (e) 

The asymptotic analysis of the solution to the second- 
order differential equation near the endpoint of the spec¬ 
trum is quite similar to that done in preceding section. 
Therefore, we only sketch its main points. 


It is convenient to define the soft-edge scaling limit as 


, 128\ 

+sy (225 j 


1/3 


Then the two-point kernel is 


1/3 


while the density of states takes the form 

1/3 

l^GUE (s) . 


V {Es) — ^N^ 


2 2 I 15 \ 

2 1-^2 I \ 


N 


32, 


(28) 


11 \ 

K {es,Es') = X%D% {—j Ague (s, s') (29a) 


(29b) 


The fitting arguments, based on the expansion of the 
Eq. (29b) and of the bulk density of states 0 


^^buik(es) - Jl 


3 














X 


3 + 4 




(30) 


near the soft edge [when 1 <C s yield = 

Dat (15/4)^^^ . Combining Eqs. (29a), (28), and (8), we 
end with the following expression for the two-point kernel 
in the soft-edge scaling limit: 


(s, s') = Kqjje (s, s') . (31) 


This formula demonstrates that in the soft-edge scaling 
limit the eigenlevel properties for the random matrix en¬ 
semble with sextic confinement potential are determined 
by the same Airy law which is inherent in CUE. 


IV. CONCLUDING REMARKS 

We have considered for the first time the correlations 
of the eigenlevels near the soft edge for two strongly non- 
Gaussian ensembles of large random matrices possessing 
unitary symmetry, and associated with quartic and sextic 
confinement potentials. Our treatment has been based on 
the analysis of the second-order differential equations for 
the corresponding “wavefunctions” near the soft edge. In 
both cases it was found that correlations between appro¬ 
priately scaled eigenvalues are universal, and character¬ 
ized by the Airy two-point kernel Eq. (3) which previ¬ 
ously has been found for CUE. 

Together with the fact of universal behavior of the den¬ 
sity of states, previously proven in Ref. |^, the consid¬ 
eration presented gives a strong impression that spectral 
correlations in invariant ensembles of large random ma¬ 
trices with rather strong and monotonous confinement 
potential are indeed universal near the soft edge. 
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